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Abstract 



A small time asymptotics of the density is established for a simplified (non- 
Gaussian, strictly hypoelliptic) second chaos process tangent to the Dudley relativistic 
diffusion. 

1 Introduction 

The problem of estimating the heat kernel, or the density of a diffusion, particularly as 
time goes to zero, has been extensively studied for a long time. Let us mention only the 
articles [V], [A], [BAl], [BA2], [L], [ERS], and the existence of other works on that subject 
by Azencott, Molchanov and Bismut, quoted in [BAl]. 

To summary roughly, a very classical question addresses the asymptotic behavior (as 
s \ 0) of the density p s (x,y) of the diffusion (x s ) solving a Stratonovich stochastic 
differential equation 



where the smooth vector fields Vj are supposed to satisfy a Hormander condition ; the 
underlying space being M. d or some <i-dimensional smooth manifold Ai. 

The elliptic case being very well understood for a long time ([V], [A]), the studies 
focussed then on the sub-elliptic case, that is to say, when the strong Hormander condition 
(that the Lie algebra generated by the fields Vi, . . . , has maximal rank everywhere) is 
fulfilled. In that case these fields generate a sub-Riemannian distance d(x, y), defined as in 
control theory, by considering only C l paths whose tangent vectors are spanned by them. 
Then the wanted asymptotic expansion tends to have the following Gaussian-like form : 




Ps 



(x, y) = s~ d / 2 exp ( - d(x, y) 2 /(2s)) £ le (x, y) s l + 0(s n+1 ) 



(1) 
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for any n G N*, with smooth 7/s and 70 > , provided (x,y) does not belong to the 
cut-locus (and uniformly within any compact set which does not intersect the cut-locus). 
See in particular ([BAl], theoreme 3.1). Note that the condition of remaining outside the 
cut-locus is here necessary, as showed in particular by [BA2]. 

The methods used to get this or a similar result have been of different nature. In [BAl], 
G. Ben Arous proceeds by expanding the flow associated to the diffusion (in this direction, 
see also [Ca]) and using a Laplace method applied to the Fourier transform of x s , then 
inverted by means of Malliavin's calculus. 

The strictly hypoelliptic case, i.e., when only the weak Hormander condition (requiring 
the use of the drift vector field Vq to recover the full tangent space) is fulfilled, remains 
much more problematic, and then rarely addressed. There is a priori no longer any reason 
that in such case the asymptotic behavior of p s {x,y) remains of the Gaussian-like type 
(1), all the less as a natural candidate for replacing the sub-Riemannian distance d(x,y) 
is missing. Indeed this already fails for the mere (however Gaussian) Langevin process 
(/3 S , Jq (3 T dr) : the missing distance is replaced by a time-dependent distance which presents 

some degeneracy in one direction, namely ^\(x — y) — § (x — y)\ 2 + ^\x — y\ 2 , see (11), 
(13) below. See also [DM] for a more involved (non-curved, strictly hypoelliptic, perturbed) 
case where Langevin-like estimates hold (without precise asymptotics), roughly having the 
following Li-Yau-like form : 

C- 1 s~ N e~ Cd ^ y)2 < p s (x, y) <C s~ N e~ c ~ l d ° {x °' y)2 , for < s < s . (2) 

In this article, an interesting case of rather natural hypoelliptic diffusion is considered 
first : that of a relativistic diffusion, first constructed in Minkowski's space (see [Du]), which 
makes sense on a generic smooth Lorentzian manifold, see [F-LJ1], [F-LJ2], [F-LJ3]. In the 
simplest case of Minkowski's space, it consists in the pair (£ s , £ s ) G M 1 ^ x M d (parametrized 
by its proper time s, and analogous to a Langevin process), where the velocity (£ s ) is a 
hyperbolic Brownian motion. In the general case this Dudley diffusion can be rolled without 
slipping from a reference tangent space to the Lorentzian manifold, see [F-LJ1]. Note that 
even in the Minkowski space, there is a curvature constraint to be taken into account, 
namely that of the mass shell H d , at the heart of this framework. Moreover the relativistic 
diffusion is never sub-elliptic, but only hypoelliptic, and a priori a Gaussian-like asymptotic 
expansion as (1) does not even make sense, since there is no longer any natural candidate 
to replace the sub-Riemannian distance d(x,y). See however [BF], where some non-trivial 
information is extracted about the relativistic diffusion, by considering the sub-Riemannian 
distance generated by all fields V , V\, . . . , V k (i.e., not only V 1: . . . , V k ). Talking of this, 
an important feature of the strictly hypoelliptic case, which is fulfilled in the relativistic 
framework, is when the graded geometry generated by the successive brackets of a given 
weight is (at least locally) constant (see [NSW], and also [T]), yielding homogeneity in the 
afore-mentioned time-dependent distance. 

To proceed, we shall compute a Fourier-Laplace transform, which seems to be the only 
way of getting any quantitative access to the density kernel of the relativistic diffusion 
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([BA1] relies already on the Fourier transform, but then the method followed by G. Ben 
Arous is based on a stochastic variation about a minimal geodesic, which does not exist 
here, and on the local strict convexity of the energy functional (due to the sub-ellipticity), 
which does not hold here) . Because of the singularity in the most natural polar coordinates, 
we first choose alternative, less intuitive but smooth coordinates, and using them, partially 
expand the relativistic diffusion to project it on the second Wiener chaos, thereby exhibiting 
a simplified "tangent process". For this simplified process the Fourier-Laplace transform is 
exactly computable. Then analyzing its inverse Fourier transform very carefully and using 
a saddle-point method allows to derive an asymptotic equivalent for the density q s of this 
tangent process, as time s goes to zero (see Theorem 6.3 below). As in the Langevin or 
in the more sophisticated case of [DM], the exponential term is given by a time- dependent 
distance, namely the same as the afore-mentioned Langevin one, the strictly second chaos 
coordinate appearing only in the off-exponent term. The initial analogous question about 
the relativistic diffusion remains open, as the degree of contact between both considered 
processes (the effective computation of the Fourier-Laplace transform being bounded to 
the second Wiener chaos) seems so far too weak, to allow to deduce a former asymptotic 
behaviour from the second one. The non-appearance of the non-Gaussian coordinates 
in the found asymptotic exponent lets however think that this could remain so for all 
higher order chaos terms of the Taylor expansion. Thence a tempting guess, resulting from 
both the sub-elliptic case (1) as solved by [BAl], the modified Li-Yau-like estimates (2) 
obtained in [DM] and the present work, whose main result is Theorem 6.3 below, is that 
an expansion having the following form could (maybe generally, under consistency of the 
Lie graded geometry) hold : 

(n x 
e=o ' 

In this first attempt we restrict to the simplest case of the five-dimensional Minkowski 
space IR 1 ' 2 x M 2 and then to its five-dimensional second chaos tangent process. The case 
of the generic Minkowski space R 1,d x M d is actually very analogous, but would mainly 
bring notational difficulties without modifying the method. We hope that this particular 
toy example will allow to understand better what can happen and could be undertaken, 
concerning small time asymptotics of the relativistic diffusion itself, and then maybe in 
some more generic strictly hypoelliptic framework. 

The content is organized as follows. 

In Section 2 are mainly described the setting and the smooth parametrization used then. In 
Section 3 the simplified "tangent process" (Y s ) to the relativistic diffusion (X s ) is exhibited. 
In Section 4 the Fourier-Laplace transform of the tangent process (Y s ) is computed. In 
Section 5 a closed integral expression for the density q s of (Y s ) is explicited. Section 6 yields 
a precise (off the diagonal) equivalent for the density q s as s — >■ 0, which is the content 
of the main result (Theorem 6.3). Section 7 contains three rather technical proofs, which 
have been postponed till there to lighten the reading. 
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2 A smooth parametrization of H 2 x R 12 = T+R 1 ' 2 

The Dudley relativistic diffusion X s = (£ s , £ s ) (see [Du]) lives in the future-directed unit 
tangent bundle TjR 1 ' 2 = H 2 x R 1 ' 2 to the Minkowski space (R 1,2 , (•, •)) (or alternatively, in 
its frame bundle, isomorphic to the Poincare isometry group V 3 = PSO(l,2) oc R 1 ' ). We 
classically identify the hyperbolic plane H 2 C R 1 ' 2 with the upper sheet of the hyperboloid 
having equation |£°| 2 — l^ 1 ) 2 — |£ 2 | 2 = (£, £) = 1 within IR 1 ' 2 (endowed with its canonical 
basis (e ,ei,e 2 )). The velocity sub-diffusion (£ s ) is a hyperbolic Brownian motion, and we 
merely have d£ s = £ s ds . The parameter s is precisely the physical proper time. 

We shall use the following smooth coordinates (A , fx , x , y , z) G R 5 on T+R 1 ' 2 : 

e° = chAch/i; ^chAsh/i; £ 2 = shA; £° = x ; ? = y ; £ 2 = z. (3) 

In these coordinates the Dudley diffusion X s = (A s , fx s , x s , y s , z s ) G IR 5 satisfies the 
following system of stochastic differential equations (for independent real Brownian motions 

w,0): 

2 dB 
d\ s = a dw s + ^- th A s ds ; dfx s = a - ^ ; (4) 

<ix s = ch A s ch fx s ds ; dy s = ch X s sh /i s ds ; cfe s = sh A s ds . (5) 
The infinitesimal generator C of (X s ) reads in these coordinates: 



^ 2 



B 2 B B 2 

+ thA^ + (chA)- 2C " 



B B B 

+ ch A ch a — — h ch A sh a — — h sh A — • (6) 

ox ay Bz 



_BX 2 BX v ' Bfi 2 

Consider the following smooth vector fields on IR 5 : 

Vl :=(T m > V * := -TTTT' V c \-= 7T + T thA m + chAch/i — + c h A sh/i — + sh A — , 
a A chA a/i as ^ a A ox By Oz 

and also V :— Vq — • We have then 

£ + |(V? + V 2 2 )+\/ ', [y 1 ,y '] = ^^ + ( jshAch / ,^ + ( 7shAsh / ,^ + ( 7chAf , 
[ y 2,Vo'] = 4^ A + ash/i^ + ach/i^, [y 2 , [V 2 , V{\] + 

and then (v^ V 1: V 2 , [V u V£[, [V 2 , V£], [V 2 , [V 2 , V$\ ) has full rank 6 at any point, so that the 
weak Hormander condition holds. 

Hence by hypoellipticity, C admits a smooth heat kernel p s (Xo; X) (s G R+ , Xq, X G R 5 ), 
with respect to the Liouville measure L , which reads L(dX) = ch A d\ dfx dx dy dz : 



E 



x [F(X s )] = / p s (X ; X) L(dX) = / p s (X ; A , fx , x , y , z) ch XdXdfxdx dy dz . 
Jt^r 1 ' 2 Jr 5 

An open question is to estimate p s (X Q ; X), for small proper times s . 
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Up to apply some element of the Poincare group V 3 , we can restrict to X = (e , 0) = 
(0, 0, 0, 0, 0). Thus we have to deal with p s (0, X) = p s {X), for X = (A , // , x , y , z) e R 5 . 

The underlying unperturbed (deterministic) process X° = (A , fx , x° , y° , z°) G M 5 
solves : 

rfA° = 4thA°rfs; /i° = 0; dx° s = ch X° s ds ; y° s =0; d z° s = sh \° s ds , 
and then is merely given by the geodesic X® = (0 , , s , , 0) (for any proper time s). 

Up to change the speed of the canonical Brownian motion (w, /3), by considering (w a 2 s , /3 CT 2 S ) 
instead of (w s , j3 s ), we can absorb the speed parameter a , and then suppose that a = 1 . 



3 A process tangent to the relativistic diffusion (X s ) 



The main Theorem 2.1 in [Ca] could apply here (beware however that Vq is unbounded), 
yielding a full general Taylor expansion for the diffusion (X s ), in terms of the above vector 
fields Vi,V2,Vo, their successive brackets, and of the iterated Stratonovich integrals with 
respect to (w,j3). Indeed, Equations (4), (5) read equivalently : 

dX s = V 1 (X S ) dw s + V 2 (X S ) d(5 s + V^X 8 ) ds . 

In [Ca], the successive remainders corresponding to the truncated Taylor expansion are 
controlled in probability. In this spirit and also almost surely, the process X s is approached 
as follows. 

Lemma 3.1 (i) For any e > 0, almost surely as proper time s \ we have : 

X s = (w 8 + \j w T dT, f3 s + o(s 3/2 - £ ), s + [ff + w 2 T ]dT, f3 T dr, J^w T dT^) +o(s 5/2 - £ ). 

(ii) Setting R' s := (x s , y s , z s ) — {^s + [(3 2 + w^dr, J (3 T dr,J w T dr^ and 
R s := (A s ,/i s ) — (w s ,j3 s ), there exist c,k>0 such that for any R > c we have: 



lim P 

s\0 



sup || i m, / 

0<t<s 



RA\>Rs 3 / 2 



< e~ RK/c and lim P 

s\0 



sup II^H > Rs 5/2 

0<t<s 



< e~ RK ' c . 



The proof is posponed to Section 7. 

Remark 3.2 More precisely, concerning the martingale (/x s ) we have 

^ = [ = [ (1 " * W * + ° {r2 ' £) ^ T =Ps ~ l2 [ W * dPr + ° (s5/2 " } - 

But the method used then does not work with the non-quadratic martingale / w 2 df3 T 

Jo 

(which equals w 2 s (5 s — 2 J* j3 T w T dw T — J* (5 T dr = o(s 3 / 2 ~ e )). 
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As a consequence, we shall use a perturbation method, approaching (for small proper 
time s) the relativistic diffusion X s by means of the M 5 -valued "tangent process" : 

Y s := (w s , fj a , \ J [(3 2 T + w 2 ] dr , J (3 T dr, J w T drj =: (w s , (3 S , A s , ( s , z s ) (7) 

which is not Gaussian, but has its third coordinate A s in the second Wiener chaos. This 
actually yields the orthogonal projection of the process (X s ) onto the second Wiener chaos. 

Remark 3.3 The fact that the second chaos term is needed in the approximation (without 
it, the tangent process would clearly not admit any density) makes a significant difference 
with the situation exhaustively investigated in [DM], where the approaching process is 
Gaussian. This can no longer be the case in the present setting, though both settings 
share the feature of being strictly hypoelliptic. A difference between both is the curvature, 
at the heart of the relativistic realm (even in the present Minkowski-Dudley flat case), due 
to the mass shell constraint on velocities. 

Note that for any fixed proper time s > we have : 

Y s l = (y/Z Wl , V~sf3i, 4 j\(3 2 T + w 2 T ]dr , V^jy r dr, ^ j\ T dr). (8) 

Denote by q s = q s (w, (3, x, (, z) the density of Y s with respect to the Lebesgue measure 
A(dY) = dwd(3dxd(dz (i.e., not the Liouville measure L) on M 2 x M* + x M 2 . 

By the scaling property (8), it must satisfy : 

q.(wJ,x,C,z) = ^ qi \^ 8 ,^,^, 7 =, 7 =), (9) 
and otherwise : q s (w, (3, x, (, z) = q s (—w, (3, x, (, —z) = q s (w, —(3, x, — £, z) . 



4 Fourier-Laplace transform of the tangent process 

4.1 Fourier-Laplace transform of the simplified process Z s 

We need information on the density at time s of the tangent process (Y s ) given in (7). 
By the independence of w,j3, it will be enough to consider the density of 

Z s :— (w s , J w T dr,J w 2 dr^ = (^/swi,Vs^J w T dr,s 2 J wldr^j. (10) 

Note that this simplified tangent process Z s does not have any component beyond the 
second chaos. Because of the scaling property (10) of Z s , its density q®{w,z,x) satisfies 

0/ \ 1 o ( w z x \ 
Qs K *> *) = ?i {-7= s ' -7=f ' 7^ J =9- *) • 



Of course, the Langevin process (w s , J w T dr^j is Gaussian with covariance 



s s 2 /2 
s 2 /2 s 3 /3 



2 m „3 /q ) ) so that it has the well-known density 



(w,z) 



71 S s 



(6z 2 ~6s zw+2s 2 w 2 )/s 3 



71 S z 



exp 



6 / s \ 2 w 2 

— [ z w ) 

s 3 V 2 ) 2s 



(in particular, the expected value of w T dr, conditionally on w s = w, equals su>/2.) 

The law of the variable Z s is not at all that simple, but it is known (see [Y], [CDJR]) that 
its Fourier-Laplace transform is computable. The following lemma is proved in Section 7. 

Lemma 4.1.1 The law of the variable Z s of (10) is given by: for any s > and real 
r, c, b, 



E 
1 



exp rw s + cj w T dr — —J w 2 dr^j 



^Mbs) 



exp 



th(frs) 
2b 



sh 2 (fe/2) bs-th(bs) 



b 2 ch(6s) 



rc 



2 6 3 



Of course, for 6 = we recover the Fourier transform of (11), namely e ( r +src+ 3 c ) s / 2 . 

Proposition 4.1.2 The x-Laplace transform of the variable Z 1 of (10) is given by: for 
any real w,z,b, 

[bz-th(b/2) w] 2 +cothb [b-2th(6/2)] w 



b exp 



e "2 x q®(w, z, x) dx = 



l-(2/6) th(6/2) 



27ry / [6-2th(6/2)]sh6 



2 7 r v /[6-2th(6/2)]sh6 



x exp 



6 2 fit; - 2^) 2 

x 



8 l-fcoth(f) 2 



^ 2 -^coth(|) W 2 



we recover 



This is of course consistent with (11), via b — > 0; and integrating with respect to dw dz , 

qi(w, z, x) dw dz dx = (ch6) _1//2 , as it must be. 

Proof We invert the Fourier transform in Lemma 4.1.1 by Plancherel's Formula : 

b 2 

e --2-x qQ( w ^ Z) x) dx 
th6 " 1 ' 2 



1 



4n 2 Vch.b Jr 2 
1 



e -V=l[wr+zc] exp 



47T 2 y/chb Jw 2 



e - V - 1[wr +zc] exp 



2b 
thb 



b 2 ch b 



2b 



r + 



th(6/2) 



c - 



2 b 3 

6-2th(6/2) 
263 



dr (ic 



— / 

s/chb Jr< 



47r 2 Vchb 



r+ (2 



th(i>/2) 



.).] 



cxp 



th& 2 &-2th(&/2) 2 



2b 



— bw z 

b e 2th& 
x — / e 



sh& 2tt 



-( 2 _i!^Z£) w ) c 



exp 



2 b 3 

6-2th(6/2) 2 ' 



2 6 3 



— 

e 2 th b 



sh& 2tt 
b 2 



x 



& 3 



&-2th(&/2) 



x exp 



b 3 



2[6-2th(6/2)] 



2tt 



27ry / [6-2th(&/2)]sh6 



x exp 



[bz - th(&/2) w] 2 + coth 6 [6 - 2 th(6/2)] 
2 [1 - (2/&)th(&/2)] 



. o 



4.2 Fourier-Laplace transform of the tangent process Y s 

We use Lemma 4.1.1 to express this Fourier-Laplace transform. 

Lemma 4.2.1 The law of the variable Y s of (7) is given by: for any s > and real 
r, e,&,7,c, 

E 



1 



ch(bs) 



cxp 



exp (y^l [r w s + q f3 s + 7 C s + c z s ] - & 2 A,) 

rc + £7) - 



tMM (r2 + , 2) _ 2 (r c + , 7) - (c 2 + 7 2 )' 
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b 2 ch(6s) 



In particular, the law of A s is given by E [exp(— & 2 A S )] = l/ch(6s). 

Proof This follows directly from Lemma 4.1.1, by independence of w and (3 . o 

Proposition 4.2.2 The x-Laplace transform of the variable Y 1 of (7) is given by: for 
any real w, (3, z, (, b, 

poo 

/ e~ b2x q 1 (w,0,x,C,z)dx 
Jo 



b A exp 



'0 



8tt 2 [6ch(|)-2sh(|)] sh(|) 



■ *«ac,*(6) • (12) 



In particular, *& w ,p£ >z (0) = ^ exp 
density of (wi , f3 1 , d , zi) . 



W 2 +f3 2 



6[z — w/2j 2 — 6(C — /3/2) 2 is the marginal 



Proof As Lemma 4.2.1 follows from Lemma 4.1.1, this follows merely from Proposition 
4.1.2 by independence of w and (5 . Indeed, for any test functions /, g on M 2 we have: 



f(w,z)g((3,() 



e 6 x qi(w, (3, x, (, z) dx 



dw dz d(3 d( 
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= E 



[ f{w,z) f 
Jr 2 Jo 



E 



e 2 x q 1 (w, z, x) dx dw dz x 



[am f 

Jr Jo 



e~ x q»(f3X,x)dxdf3d( 



p r /»oo 2 /*°° 2 

= / f(w,z)g(p,{) / e-^ x q°(w,z,x)dxx e~^ x q°{(3, (, x) dx 
Jr 4 Uo Jo 

and the claim follows directly from Proposition 4.1.2. o 



dw dz d(3 d( , 



Lemma 4.2.3 All solutions z e C* of the equation z = thz belong to the imaginary 
axis, and form a sequence 1Z = T y n \n E N}, with ^ < y < y 1 < . . . < y n f - oo . 

Proof Writing z = x + T y £ C*, we have z = thz e 2z = and then equivalently 
cos(2y) = ( ]zf-£ 2 e~ 2x and sin(2y) - 



-2x 



whence 



e 4* ((1 _ a;)* + ^ = 4y 2 + (1 - X 2 - y 2 ) 2 . 

The latter is equivalent either to x = , or to y 2 = — (1 — x) 2 (which is excluded), or to 

Then using this last value of y 2 , by the above we must also have 



< 2 < 



sh (2a) 



sin(2y) = - sh (2x), which is impossible since for any i,i/6R* we have 
and clearly z cannot be real. Hence we are left with z = ±y/—l y , with y > and then 
y — tgy . The claim follows, with moreover ((n + 3/2) 7r — y n ) \ , and ^ > y > ~nf 
since tg 



17n 
12 



cotg ^ = 2 + v / 3<^- o 
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The Laplace transform ^^^(fr) in Proposition 4.2.2 is a meromorphic function of b G 
C , with singularities at the points of 2rr Z* and at the non-null zeros of [b— 2 th(6/2)] , 
that is to say at the points of 271 (according to Lemma 4.2.3). 



4.3 Complement: Density a s (x) of the variable A s 

Denote by a s = a s (x) the density of the variable A s , so that for any s > we have 
a s (x) = J q s (w,f3,x,(,z)dwd/3d(dz = ^ai(^), by (9). 

Lemma 4.3.1 The density a\ is smooth and bounded (with bounded derivatives), and we have 



4 r 

(*) = " / 



cos (2x y 2 - y) 4 f°° cosy cos (2xy 2 ) 

—To — ; o—yshydy+- — ■ = — 

sh z y + cos^ y it J sh z y + cos^ y 



e y ydy . 
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Proof According to Lemma 4.2.1, for any 77 > we have: a±(x) = 

dy = — / — ^=^== - dy . 



2vr . / _ 00 



Expanding and letting rj \ we obtain : 

. A f°° chy cosy cos(2xy 2 ) + shy sin y sin(2xy 2 ) 

a 1 {x)=- -3 — ■ 5 ^ 

vr y sn 2/ + cos^ y 

4 Z* 00 cos(2xy 2 - y) 4 f°° cosy cos(2xy 2 ) _„ , 

= _ / — ui — : — 2 — v sh -ydy + - —-k — ■ 5 — e v d v- 

7r Jo sh^y + cos^ y 7r y sh z y + cos^ y 

4 Z" 00 ychydy 
As a consequence, ai is smooth, and bounded by — 



7r y o sh 2 y + cos 2 y 

Remark 4.3.2 We have sh 2 z + cos 2 2: = O ch 2 z = sin 2 z O z G (1 ± v/-l)^(l + 2Z), 

i.e., |z| G 7 Argz = j modulo | • 

Proof There is clearly no real solution, since chx > 1 > sinx for any x£R*. Consider then 
z = (\/— T + a) x , with x,a G R. Up to change z into — z, we only have to consider the 
equation ch z = sin z . Then 

chz = sinz -4=>- ch(ax) cosx = chx sin(ax) and sh(ax) sinx = shx cos(ax) 

=>■ ch 2 (ax) cos 2 x + sh 2 (ax) sin 2 x = ch 2 xsh 2 x -4=> (ch 2 x — cos 2 x) (sh 2 (ax) — sh 2 x) = . 

As x = is not a solution, the only possibility is sh(ax) = ±shx , whence a = ±1, and then 
cos x = a sin x . This yields the claim, o 

( 1 "\ 2 2 

Proposition 4.3.3 We have ai (x) = 2vr (-l) n (n + \) e'^) w x , so that 

neN 

7re~ 7r2x / 4 (l-3e- 2w2x ) < «i(x) < Tre" 7 ^ 4 for any x > ir~ 2 , and then 
Qi(x) = 7T e~ n X I A (l — 0{e~ 2n x )) as x — > oo . Moreover, a.\ decreases on [3"7r~ 2 ,oo[. 



Proof This results from Lemma 4.2.1 and ([BPY], Table 1 continued and Section 3.3, in partic- 
ular Formula (3.11), with C\ = 2A\, got merely by expanding l/(cht)). For any 7r 2 x > 1, this 

alternate series has decreasing generic term (n + rj) e~( n+ 2") n x , whence the estimate. Finally, 

3 ( i\ 2 2 

the same holds for the series (n + i) e-y n+ 2> * x yielding a'^x), thereby garanteing a'^x) < 0, 
as soon as 7r 2 x > 3. o 



Lemma 4.2.1 and Proposition 4.3.3 at once entail the following. 
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Corollary 4.3.4 A 



„A x 



a\{x)dx defines an analytic function on 3i(A) < 7r 2 /4, equal 



^7= for < A < V Moreover, we fcove a a (x) = 4e"" 2 " /(4s2) (l - 0(e~ 2 ^ 2 x/s2 )) (as 

'S v A s \ / 



to 

f -> oo). 

Proof By Lemma 4.2.1 the set of b £ C such that / 

Jo 



_, 2 1 
e x cti(x) dx = — — - contains R. By 
v ' ch6 



/•oo 2 I 
Proposition 4.3.3 b h-> / e _b :E ai(x)(ix is analytic on < tt/2}, and so is b i->- too. 

/•oo ^ 

Hence we have / e b x oii(x)dx = for |§R(6)| <tt/2. o 

Jo cosb 

Remark 4.3.5 1) The law of 2 is that of inf {r > | \j3 T \ = l} and also that of 

^max{|/3 T | | < r < l}^ (See [BPY], Table 2, which also exhibits two other random variables 
having the same law as Ai). 2) Note that for 3?(A) < 7r 2 /4 (letting sj~- > on R+) : 

/*00 p p 

J J e Xx qi(w,/3,x,C,z)dwd(3d(dzdx = J ^ Wj ^ z (V—\) dw dfi d( dz . 



cos y/X 

5 Integral expression of the density q s (w, x, £, z) 

Let 



q s = q a (w,p,C,z) 



7rV 



cxp 



3[(sw-2z) 2 + (sf3-2() 2 } w 2 + p 



2 s 3 



2s 



(13) 



denote the marginal Gaussian density of Y s := yw s , (3 S , J /3 T dT,J w T dTj (according 
to (11) or Proposition 4.2.2)). 

Notation Consider the function $(A) = $«,,^,^, Z (A) := tyw,p,c,z(V— A), derived from 
Proposition 4.2.2. We systematically use the usual determination of the complex square 
root, cutting C along the negative real semi-axis and letting \f~- > on W* + . By the 
expression (12), we have: 



(A) 



A 2 exp 



8tt 2 [2sin(^) - v / Acos(^)] sin(^) 
Lemma 5.1 The function A h-> <&u, J( a,e, z (A) is analytic for 5ft(A) < 47T 2 . 



(14) 



Proof Note that the functions sin(^) / ^ and cos(^) are plainly analytically contin- 
ued for any A G C, and that by the expression (14), $(A) is analytically continued at A = 0, 



and analytic at any A G C* such that sin(^) ^ and tg(^p) ^ ^ , hence, according to 
Lemma 4.2.3, at those A G C not belonging to the sequence {47r 2 (n + l) 2 ,4y 2 | n G N} C 
[47r 2 , oo[ . This shows the analyticity for 3?(A) < 47r 2 . o 

The following is proved in Section 7. 



11 



Proposition 5.2 For any (w,/3,(, z ) we have 

PCX) 

$™AC, z (A) = / e Xx qi (w,p,x,(,z)dx , for 3?(A) < 4vr 2 , 
./o 



and 



e rx qi (w,(3,x,(,z)= e -^ tx <5>{r + V^lt) — , for x > and r < An 2 . (15) 

J — CO 



271 



By scaling and using (15), for any s , x > , r < An 2 and (w, /3, (, z) G K 4 we have : 

(a a \ 1 ( W X C * 



-rx/s 2 roc 



2tt s 6 



( w 13 C z "\ ' 

2 



r + V-l i/) oh/. 



Taking merely r = 0, and setting := e vizy/s $ 



for convenience, for any s,x > and (w,/3,£,z) G M 4 we have: 



1 /"°° ~ 1 
q,(w,P,x,C,z) = 7 ^— R \ * 8 (y)dy = 7 ^ / (V s (y) + V s (-y)) dy . (16) 



2tt s 6 



Now according to (14) we have : <3> 



2tt s 6 



«L,^_,_C_,^L_ 



2 3'. y s-' 2 - v cotg ( y ) ( W 2 +/3 2 ) 



8tt 2 [1 - cos( v / v /Z Ty) - ( v / 7^/2) sin( v / 7 r T^)] 

in which we have set 
5 



cxp 



B 2 V^iys-' 



l - 



tg(^) 



2 _ (gw - 2z) 2 + (gg - 2Q 2 s , = 4 W z + 4/3C-s( W 2 + /3 2 ) ^ 



Let us now write out a more tractable expression of ty s introduced above. 
First, for any real y, 9 , setting £ := -y/^-) we successively have: 

2 



vV^Ty = (l + sgnfc/JV^T) £ 



(l-sgn^v^T)/^; 



tg((i + ✓ z T)9) = ^'t^y ; + x/-T)0) = ^ - ^ sh < 29 » 



ch(20) + cos(26>) 

+ f y/V=Ty \ _ sing + sgnM-y/^Tshg f y/V=Ty \ _ sing - sgn^y^T shg 

tg l^^ ) ~ chg + cosg 5 cot S(, 2 j - — - — - 



ch(20) - cos(26>) 
ch g — cos £ 



12 



y^T^N _ (sh£ + sin£) + sgn(j/)^/^T (sh£ - sin£) 



) 



: ) ^(che + coso 

£ (ch g + cos g) - (sh £ + sin £) + sgn(y) ^/^T (sh g - sin Q 
i (ch £ + cos - 2(sh i + sin £) + 2^ (ch £ - cos 



c _ shg+sinA 2 . / shg-sing x 2 
^ ch£+cos§ J \ ch£+cos( 



=: U(y) ; 



V / \ /r T jj cotg f ^~^y \ = (shg + sing) ~ sgn(y)V^l (shg - sin£) = _ y ^ 



2 ; 2r 1 (che-cos0 
cos^a/\7^T y ^ = ch^ cos£ — sgn(y)\^-T sh^ sin£ 



sin 



— 1 y J = ch £ sin £ + sgn(y) y — 1 sh £ cos £ ; 
l-cos(V7=Ty) - (y/y/^ly/2) sinfy^y) = 
1 — ch£ cos£ — f (ch£ sin£ — sh£ cos£) + sgn(y)V _ 1 sh£ sin£ — §(ch£ sin£ + sh£ cos £) 



Set 
with 



1 - cos {VV^ly) - ( v / 7 r Ty/2) sin (>A/=Iy) 



and 



Hence 



^(0 := 



ch £ cos £ + | (ch £ sin £ — sh £ cos £) — 1 



(ch£ - cosO [(chf - cosO - £ (sh£ + sinO + ^(ch£ + cosf)] 

sh £ sin £ — | (ch £ sin £ + sh £ cos £ ) 
(ch£ - cos£) [(ch£ - cos£) - £ (sh£ + sin£) + f (ch£ + cos£)] 



V 2 F{y) 
8tt 2 



exp 



B' s - x w 2 + (3 2 



2s 



V(y) + ^ t%) 



(18) 



8tt 2 



exp 



B 



B' - x w 2 +p 2 



V=I Hf + -h- + K(0 v - ^ ^(0 



w 2 +l3 2 

s 2 ' 2s " JVS 7 y s 3 ^^""27" 



K(0 



with 



t ( t _ shg+sinA 
^ \^ ch§+cos£ J 



C sh£+sin£ 
^ ch£+cos£ 



+ 



sh£— sin£ 
ch£+cos§ 



2 ' 



[/,(£) := 



o t3 ( shg-sing \ 

l^ch^+cos^y 



sh§+sin£ 
ch£+cos§ 



ch £+cos £ 
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4f (chf -cosf) 
Using these auxiliary functions, (16) reads 



2(chf-cos0 ' 



q s (w,p,x,C,z) 



qs{w,i3,c,z) 

(3tt/2) s 2 



i r 

8tt^ J 



3? 



F(y)e 



B' s - X , W 2 +fi 



3? 



F(y)e 
3n s 2 



2C L 



e -$Vr(t)-^v r{0y 2 dy 



e -|f[l/ r (0-12]-!£^[V r (0-l] 



qc/ r (0-i2] 



[Vr(S)-l] t5 



(19) 



with 



= 2F r (0 cos(A s (0) - 2F,(0 sin(A s (0) , 

where we set A s (0 := 2£ 2 (§ f/,(0 + ^ + ^ ^(0)- 
Note that the functions F r , F i: U r , U i: V r , V i: G s are all even. 



6 Small time asymptotics for the density q s (w, (3,x,(, z) 

We shall here use the expression (19) for q s (w, /3, x, (, z), to derive its asymptotics as 
s \ 0, proceeding by adapting the saddle-point method (see [Co] for example). This will 
require the following asymptotics for the auxiliary function entering that expression. 



6.1 Auxiliary asymptotics 

As ^oo wehave U r (£) = 2£ + 4 + | + 0(±) ; U^) = 1 + 0(\) ; K(0 = |[i+o(e-«)] 

vi(0 = 1± ^ =!l ; iW) = o(*±). 

Then near we successively have : 

(S) 2 +(f(i-l« 4 )) 2 +f(n S(i-^)+o«") 35 

rr (£\ 30 (-*• ~ 648^ + ^£ )) 6 _ 2119 £-4 , /-/Vt8\ . 

I(i -i^+w)) 5 18900 4 lu ' 



14 



f* - Hi + (^ 4 ) ; = ^ + o(^) ; 



a.(0 



3 



= 2/1. f- 
where we have set 



2119 i 1 w 2 +f3 2 
9450 IF "I" 63 12 s 



[e 6 + ^ 10 )] = 2^ 2 [l + 0(£ 4 )], 



/is := 



5 s 3 



+ 



B' s — x w 2 + (5 7 



+ 



24 s 



Let also za := 



17 R 2 w 2 + (5 2 



35 s 3 



+ 



360 s 



(20) 



so that 



e -^[^)-i2]-^h±£_[K(0-i] 



s? ri7 



[ilo« 4 + « 8 )] = e-M^+OK 8 )]. 



It is not difficult to see that actually > 1 and U r (£) > 12 , for any real £ . 

Furthermore, 



tf r (0-2e-4 = 



4£(1- 



2(e £-sing) \ _ g / chg-cosg \ _ r> f 3 / e^+cos g+sing \ _ ^ f 2 / e^+cosg-sing 
ch£+cos£ / "Vchf+cos^ Z ^ I ch£+cos£ ) ^ ? V 



ch£+cos£ 



ch £+cos t 



) 



(p — 1^2 _i_ 1 I rW e-€+cosg-sing \ _ 4 cosg 
^ T J- T ^ ^ chg+cosg ^ ch£+cos§ 



> 



4£-8- 



4(g+l) 3 
ch£+cos£ 



(e-l) 2 + l + ch?+cosf 

for £ > 27T, and similarly 



4g-8+Q(g 3 e-£) 
£ 2 -2£+2+0(£ e -e) 



^ 4£-8-9(£ + l) 3 e-g ^ 



4£-9 



2 2tt-9 
> - x 



(£-l) 2 + l + 9(£ + l)e-« £2_ 2 £ + 2 .2-£ 7T-1 + 



li 

20?r 



rlU 2 2(chf-cosf) ch£-cos£ 



s 20 



so that K(0 > vr - ^ > 3 for £ > 2vr. 

Therefore, in a small neighbourhood [27r — e, oo[+v^-T [— s,e] of [2ir, oof we shall have 

»[t/" r (0 - 2f] > 4 and 9fc[V r (0] > 3. (21) 
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6.2 Changes of contour and saddle-point method 

Note that by Lemmas 4.2.3, 5.1, (14) and the changes of variable: A = \J —\ y = 
±2\/— 1 £ 2 , the poles of the integrand in (19) are located at e*^ =Tfc7r//2 (l + \J — 1 ) y n and 
e v/ ^ Ifc7r / 2 (l + v / — T)(n + 1) 7T, with n, k G N and ^ < y < y 1 < . . . , so that this integrand 
is analytic at with convergence radius y/2ir. 

In particular, we may perform the following changes of contour in (19) : 

.2 q s (w,(3,x,(,z) _ r nf n-4luA0-i2]-^\y r ®-i]t5. 



q s (w,/3,(,z) Jo 

= f^^^ " ( Fr(0 + F . (0 ) e v^T A s(0 -f [^-121-4^(0-1] ^ ^ 

Jo 



+ / (F r (0 + V^T F,(0) 6^^(0-^^(0-121-^^(0-1] ^5 ^ 

+ / (F r (0 + v^T F,(0) e^ 1 MO-^O-iq-^WO-i] ^ ^ 

V(1+V=T) M - "+2tt 

H 1- ' V^-T)/u7' , +0O , B 2 ,„2 + «2 

+ / ^ (F r (0 - v 7 ^ m)) e'^ 1 A ^-^ u ^-^-^^-^ 



'(l-/=l)/i s "+27T 

=: ^ + ^? + J7 + ^ + jr' + ^ , 
where 77 > | will be specified further. 

Note that the estimate (21) and the control F(2£ 2 ) = O(^) ensure the vanishing of the 
unmentioned limiting contribution (for any large enough fixed s, provided lim fi s = 00) in 

s^>0 

the above changes of contour : 
lim 



/\K±(1+V--UM« , R 2 ,„2_l_«2 

im / (F r (0 + v^F t (0)e^ As «)-^^- 12 ]-T[^M] e 5^ = Q 

^°°Jr 

Now on the one hand, setting £ = (1 + ^/— T)^^ t by the above we have: 
Jo 
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5/x 



6t; 



./O 



1 + 



15^/^1 i 620659 V /3 T 



+ 



2/ i - 4 "^ ' 27120 



+ 0(4 



-4jj\ 



di 



(i + (/^ + ^(/O)) 



" 5/4" 
so that for \ < rj < \ : 



1 + + + 0(^)\ t 5 dt , 



jO jo = J! 



(i + e>(^: 4r 0) [\-^~ 2Vt2 t 5 dt 

Jo 

, r 1 \ „.. 2J .. i + 0(^- 4 ") 

io 



e "u dti = 



20 ^ 



1 + 0(^7^ 

provided lim /x s = oo and z/ s = 0(fj, s ). 

Then we use the saddle-point method (as described in [Co]) to deal with the interme- 
diate part (JJ + JJ). By the above we know that we have 

u t (0 = 1 + ^ m) , vs(0 = h + e m) , u r (0 - 12 = £ 4 u r (0 , v r (0 -1=^ v r (0 , 

with even functions Ui ,Vi ,U r ,V r that are analytic outside the sequence 

jgv^T _|_ ^JZZ\ ) y n | n, A; G N, ^ < y < Hi < ■ ■ ■ } and then converge in the compact 

disc (centred at 0) of radius > 2tt. (Note that the proof of Lemma 5.1 shows that the 

points e^^/^l + v^X™ + 1) 

7T are poles only for the functions F T) F { .) Therefore we 

can write the phase as follows : 



Setting M : = max {HOI | £ e C, |f | = by Cauchy's inequality, for |^| < vr + ^ 

we have 



-i M0 - ^[^(O - 12] - ^[K(0 - 1] 



k(0-2v^T/i s e 2 | < ^ 



' 289tt 2 



< 



4(72) 2 M|£| 



|£| 2 ) ~ 289(579 -408^)^ 



< 1156 M. 



72 V 72 

As before near we have (-F r (0 + v^— T -Fj(0)£ 4 = 6 + 0(£ 4 )- Hence for small s 



J s 



f [e + ©(Ki+v^K" + t| 4 )l e 2 ^ T ^(( 1 -^ T ^' ,+ *) 2+0 ( 1 )((i+v / ^T)/i;' ? + t) 

Jo L J 

= 0(1) r e'^-^dt = 0{e~^ l^p). 
Jo 



dt 
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The same of course holds for JJ. 

To deal with the remaining contribution (J s °° + J^°), we use the lower estimates near 
infinity: (21) computed above, as follows: 



(l + V^T)Hs V +OD 2 2,^2 

(i? r (0 + v/^T F,(0) e^ 1 ^5 di 



roc 
J2-K 



o 



t 



e <>■ 



slot fi i «i 2 +/3 2 r 8fl2 «, 2 +/3 2 f 00 r 2jB 2 n 

J2n 



0(1) e 



8B 2 «) 2 +/3 2 2B 



+ 1) 



2tt 



e-* t 4 c/t = 



-4(tt-2) 



Bj w 2 +/3 2 



and the same of course holds for J°°. 



So far, we have obtained : 



^ s2 q s (w,(3,x,(,z) = J + J0 + j; + J; + J oo + joo 



g s (w,/3,C,^) 
1 + 0(u s ^) + ^-4^/^ + 0[(^i + 

l + 0(/iJ- 4 ") „ ! 1 

- for i < 77 < | , 



20 a*: 



-4 (7r -2)$-^ 



20 ^ 

provided both lim fi s = oo and z/ s = 0(^ s ). 

By (17) and (20), this condition holds as soon as both lim jj, s = oo and < 2 + £ +— (for 
some e > 0), and then also as soon as both \imu s = oo and 2s x < (z 2 + ( 2 ) +es 2 (w 2 + (3 2 ). 

Finally, under this condition, as s \ and for any positive £ , we have obtained : 

1 + 0(/i e_1 ) 

q a (w,p,x,C,z) = —— r^g- x q a (w,p,(;,z) . 

60 n s 2 n 6 s 

The result of this section (and main result) is thus the following off-diagonal equivalent. 



Theorem 6.3 As s \ 0, for any positive e, uniformly for x > and (w, (3, (, z) G 
such that: 



3[(^-H) 2 + (C-if) 2 ] , w 2 + P 2 x x z 2 + ( 2 u s 
l^s = - o 1 — ~ — > oo and — < 1 , (22) 



5 s 3 



16 s 



s 2 s 3 



we have 



1 + Oinir 1 ) 

q 8 (w,0,x,C,z) = s x exp 

20 n 3 s b [i 3 



3 [(sw - 2z) 2 + (s/3 - 2C) 2 ] w 2 + P 21 



2 s 3 



2s 
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An alternative condition (to (22) above) guaranteeing this asymptotic equivalent is 



{sw-2zf + {sB-2Q 2 w 2 + B 2 , . 2 a2n 2 . 2 n2N 

7 v — + >oo and 2s x < {z 2 + ( 2 ) + es 2 {w 2 + B 2 ) . (23) 



Remark 6.4 Theorem 6.3 rather precisely yields the small time asymptotic behaviour 
of the heat kernel of the second chaos approximation (Y s ) to the Dudley relativistic 
diffusion (X s ). But this is not enough to derive any small time asymptotics for the 
density p s (\, fi,x,y, z) of the original process X s = (A s , /i s , x s , y s , z s ), even for fixed 
(A, /i, x, y, z). The reason is that the computations performed in Section 4 above cannot 
work beyond the second chaos, so that Section 3 cannot yield a sufficiently precise control 
on the gap between both tangent processes {X s ) and (Y s ). To be more specific, Section 3 
and Theorem 6.3 heuristically yield : 

p s (w,0,x,C,z) « P[A S = w, jj, s = /3, x s = x, y s = C, z s = z] w 



P 



K , B s , A s , (s, zs) = (w, B, x, C, z) + 0(R)(s 3/ \ s 3 / 2 , s 5 / 2 , s 5 ' 2 , s 5 / 2 )] + 0(e~ RK l c ) 



exp 



',[(sw-2z+0(s 5 / 2 )) 2 + (sl3~2(;+0(s 5 / 2 ))' 



2s s 



2s 



20 7T 3 s 6 n, (w +C(s 3 / 2 ), 8 +C(s 3 / 2 ), x +C(s 5 / 2 ), C +C(s 5 / 2 ), ^ +C(s 5 /2)) 

\2 

2> 



= exp 



3[(sw-2z) 2 + (sB-2Cf] _ w 2 + B 2 _ 1/2 
2*3 2s +UlS j 



+ 0(e- RK / c ), 



which were not too bad only if the additive term 0(e RK / C ) were not there to ruin such 
estimate. Indeed, even taking R x s~ 7 would only control this correction term by (at best) 
e~ s 27/3 , which would be significant only for 7 > 9/2, so that the remaining information 
would then reduce to nothing. 



7 Proofs of some technical results 

We gather here the rather technical proofs of Lemmas 3.1 and 4.1.1 and Proposition 5.2. 
Proof of Lemma 3.1 (i) Equation (4) entails that for small proper time s we have : 

\ s = w s + \( th.{w T + o(r))dr = w s + \ f w T dr + o(s 2 ) = w s + o(s 3/2 ~ £ ^ 
Jo Jo 



= w s + 
Then 



\f th(w T +o(T^ 2 - £ ))dr = w s +\ [ (w T +o(r 3 / 2 - £ ))dr = w s + \ [ w T dT+o(s^ 2 - £ ). 
Jo Jo Jo 

chA s = ch[w s + o(s^ 2 ~ £ )] =l + \ w 2 s + (s 2 - £ ) = 1 + o(s 1 - £ ) , 
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and by Equation (5) we have : 

z s = sh\ s = sh[w s + o(s 3/2 ~ £ )] =w s + o(s 3/2 ~ £ ) . 



Similarly, 



1*8 



(Also, z s = J* a/1 + i 2 dw T + J* z T dr 



= r chK^)] ^-'^ = ft + °< 

ch fi s = 1 + o(s 1 ~ £ ) , sh fi s = fi s + o(s 3/2 ~ e ) . 
(Also, sh » s = j; dP T + /„• ^ dr .) 

Finally, the result follows at once from : 

x s = f [1 + \ w 2 T + o(r 2 - £ )] [1 + \ f3 2 T + o(r 2 - £ )] dr = s + \ f [/3 2 + W 2 ] rfr + o( S 3 ~ £ ) ; 
Jo Jo 

y,= f [l + o(r 1 - £ )}[P T + o(T 3 / 2 - s )}dr= [ S p T dr + o(s 5 / 2 - £ ); z s = f w T dr + o{s^) . 
Jo Jo Jo 

(More precisely, z s = w T dr + Jq JqW u dudr + \ J^w 2 dw u dr + o(s 7 ^ 2 ~ £ ) .) 



(ii) For any s G [0, 1] we almost surely have : |A S | < sup \w\ + | / |A r | dr , whence 

[0,s] Jo 

by Gronwall's Lemma : sup |A| < sup \w\ x e s/2 . Then for < s < 2 log 2 : 

[0,s] [0,s] 

A s :=sup|A — w\ <|/ sup \w\ x e T ^ 2 dr < / sup|w|rfr. 
[o,s] Jo [o,t] Jo [o,t] 

Hence for i? > 1 and < s < 2 log 2 : 



P 



sup |A t | > i?s2 

0<t<s 



Then f3 s -pi s = 2 



P 



< P 



s s h 2 (A r /2) 
ch A T 



sup \w\dr > R 

o [0,r] 



sup |u>| > i? 

[0,1] 
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i i 3 

sup m t — fiA > Rs? 


< 2E 


_0<t<s 





< p 

s sh 4 (A T /2) 
ch 2 A T 



< 4P[wi > R] < 2e~—. 



dr 



exp - 



, so that 

s sh 4 (A r /2) 
ch 2 A T 



< 2E 



2„3 



exp - 2iTs 



A 4 dr 



< 2E 



exp — 



R 2 



sup | it; I 

[0,1] 
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2 / P 



R 2 



8 sup \w\ 

[0,1] 



> logy 



dy = 2 P 



sup \w\ > 

[0,1] 



-R 2 



dy 



< 4 



f 1 _i f—gr f 1 -r r c 

/ e 2V 81o g y dy = 4 / e V2 1og(l/y) g^/ = 4 / 

Jo io io 



Then 







z s — w T dr 




Jo 


Jo 



8 logy. 

e^'dt < 132e- R2/a /( 32 ) 1/3 . 



= jf shw r -w T + 2sh[^f^] ch[^±^] dT<jT (i? 3 r 3 / 2 + 2 sh [^-] ch[2i?v^])rfr 
< [ (R 3 + 2R) r 3j2 dr < R 3 s 5/2 for < s < s R > , with probability 1 - 0(e~ R2/2 ). 



Similarly 



/ 

Jo 



P T dT 



dr 



[ ch(w T + (A r - w T )) sh(/3 T + (/i T - f3 T )) - (5. 
Jo 

= [ (l + 2sh 2 [^] + 2sh[^f^] s hp^])(sh/3 T + 2shp^] ch[^±^]) -(5. 
Jo ' 

< jT (r 3 t 3 I 2 + Rr 3 ' 2 + OiRr^dr < R 3 s 5 ' 2 for < s < s' R > , 
with probability 1 - ^(e^ 273 / 4 ) . Finally in the same way we obtain : 

i+\J\p 2 T +wl]dr) < ^ S |ch( WT + (A r -^ r ))ch(/3 r + (/i T -/3 r ))-l-|[/3 T 2 + W 2 ] 

(chuv + 2sh[^f^] sh[^±^])(ch/3 T + 2sh[^^] sh[^±^]) -l-±[[3 2 + w 2 ] 
< [ S (|2sh 2 [^] + |2sh 2 [-] - f | + 0(i?V))dr 



dr 



dr 



dr 



< J [R a t 2 + 0(i? 2 r 2 ) Jdr < R 4 s 3 < s"' 2 for < s < s" R > , 
with probability 1 — 0(e _R2/3 / 4 ) . In particular we can take k = 2/9 in the statement, o 

Proof of Lemma 4.1.1 Let us use ([Y], Chapter (2) "The laws of some quadratic func- 
tional of Brownian motion"), considering for any b > the exponential martingale 

Ml : = exp ^ — | [w 2 — w 2 — s) — y J w 2 dr j = exp ^ — b J w T dw T — y J w 2 dr j , 

and the new probability ¥ b having on J 7 , density with respect to P. As noticed in ([Y], 
(2.1.1)), Girsanov's Theorem yields a real (P 6 , T s ) -Brownian motion (B u , < u < s) such 
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that w u = Wq + B. 



u -b w T dr , 
Jo 



which means that under P b , w has become an Ornstein- 



Uhlenbeck process, alternatively expressed by w u = e 



—bu 



( tt '° + I 



e bT dB T . Therefore 



En 



exp ( a/—T / (a + c r) dw T — y / w 2 dr 



= e - bs/2 xE b exp(^/-l J (a + cr)dw T + bw 2 s /2 
on the one hand, and on the other hand for any test-function / on IE 

dB ll -be- bu ( [ e bT dB T \du 



\ f(u)dw u = / f(u) 
Jo Jo 

= f (/(r) - be bT £ f(u) e- bu dv^dB T , 



so that 



En 



/(t) dw T 



Taking /(r) = a + cr , we obtain 



E 



aw* + c rdw 



f(r)-be br J^ f(u)e- bu duj 2 dr. 



a + cr+(a + c/b) (e b(T - s) -l) + c(s e b(r - s) - r) 



dr 



',-2 



(ab + c + bcs) e b{T - s) -c 



dr 



= b~ 2 [ \(ab + c + bc s) 2 e 2b{T ~ s) -2c(ab + c + be s) e b{r ~ s) + c 2 ] dr 
Jo L - 1 

~ 3 \{ba + (bs + l)c) 2 (1 - e~ 2fes ) - 2c (6a + {bs + l)c) (1 - e~ bs ) + fee 2 ] 

l- e - 2hs 2 , bs-l+2e- bs ~(bs+l)e- 2bs , b 2 s 2 -3+4 (bs+1) e- bs -(bs+l) 2 e - 2bs 2 



2b 



b 2 



ac + 



2fc 3 



Th 1S yi elds the covanance matrix of the ^-Gaussian variable (»,. hence its 
joint law. Namely the covanance matnx under K of (V2b Ws , VW^rd^) is 



1 - e- 2bs bs-l + 2e- bs - (bs + l) e - 2fes \ 

bs - 1 + 2e- 6s - (bs + l)e- 26s 6 2 s 2 - 3 + 4 (6s + 1) e - bs - (bs + l) 2 e - 2bs ) ' 
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and its determinant is 5 bs '■= 2(bs — 2) + 8e bs — 2(bs + 2) e 2bs , which increases with 
bs > and then does not vanish. Therefore the density of (^/2bw s , V2tft J r dw^j is 



(u,v) 



exp 



bs 



-1 



2 5, 



bs 



(a bs u 2 + 2 lbs uv+(l- e- 2bs ) v 2 } 



with a x := x 2 — 3 + 4 (x + 1) e x — (x + l) 2 e 2:E and 7^ := 1 — rr — 2e x + (x + l)e , 
so that 5 X = (1 — e~ 2x ) a x — 7 2 . Hence 



1 



27rV^ Jk 2 



En 



exp^-\/^T y (a + cr) duv + bw 2 /2 



exp 



2 5, 



(« bsM 2 + 2 7bs Mt;+ (l-e- 26s ) v 2 ) 



du dv 



^2tt(1 - e- 2fes ) 



£VEL n J ^ „ f cy/=T 7ft,, ~\ 2 , / l.^iM.") 



g4(l-e- 2i »)6 3 

v/2tt(1 - e- 2bs ) 



( ay^T 



ey^T 7b 3 
( 1 _ e -26s) v ^3 



du 



-Jbs £ 



^2tt(1 - e" 2bs ) 



"63 ' 



2 e 4(l- e -26«) 6 3 -1 1-e- 2 ^ f ■ c 7 ft V 



v /2(i + e- 2fes ) 



A/ch(fc 



exp 



3 bs/2 



A/ch(fc 



exp 



46 3 (1 + e 



-1/ S bs C 2 n l -Z- 2bS ( u_ JYbsC_\ 2 \ 

4& 3 V 1 - e~ 2bs 1 + e~ 2bs V 1 - e~ 2bs J J 
l —— (2 (1 - e- 2bs ) b 2 a 2 + (2a hs - 5 bs ) c 2 - 4 6 76s oc) 



Therefore for any b,s > and real a, c we obtain : 



exp^v^— T J (a + CT)dw T — Y J w 2 dr 



exp 



-1 



46 3 (1 + e- 2bs ) 



(2 (1 - e- 2bs ) b 2 a 2 + (2a bs - 5 bs ) c 2 - 4 lbs b ac) 



Then taking a = r + c s , we get : 



En 



exp v^T 



ru? s 



+ c / u> T dr — y w 2 dr 
Jo J ./o 
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cxp 



exp 



exp ^v^T J (r + cs — ct) dw T — y J w 2 dr^j 

2(1 - e' 2bs ) b 2 (r + c s) 2 + (2a bs - 5 bs ) c 2 + 4 lbs b(r + c s) c 



Ab 3 {l + e- 2bs ) 

(1 - e- 2bs ) b 2 (r 2 + 2s rc) + (bs (1 + e- 2bs ) - (1 - e - 2bs ))c 2 + 2fe 7bs rc 



2 6 3 (1 + e- 



2bs^ 



exp 



Mbs) r 2 _ bs - th(fe) ^ _ 2 _ 



sh 2 (6s/2) 
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2 6 3 



b 2 ch(6s) 



. o 



Proof of Proposition 5.2 Denote by An G M+ the abscissa of convergence of the 

integral, so that the map A i— >■ / e A:c gi(u>, /3, x, £, z) dx is analytic on (3?(A) < An}- By 

Jo 

Proposition 4.2.2 and Lemma 5.1 it is equal to $(A) for 3?(A) < min{47r 2 , Ao}. Hence, for 
any real A < min{47r 2 , Ao} and t G H. we have 

POO 

$(A + v /Z Tt)= / e ^ tx e Xx q 1 (w,0,x,C,z)dx. 
Jo 

Let us show now that t h- >■ $(A + V — 1 belongs to L 1 fl L 2 (R), in order to inverse the 
above Fourier transform. Of course we have to deal here with the large values of \t\ , i.e., 
of | A + v^—T t\ in the expression (14) : 



A 2 e B'\-^cotg(^)(wl+f3 2 ) 



8tt 2 [1 - cos(v / A ) - ( VA/2) sin(v / A )] 



cxp 



B 2 \ 



in which A is to be replaced by A + \/— T t, and we have set B 2 := ( w ~ 2z ) ±<^jC) an( j 

v 4wz+4/3C-m> 2 -,8 2 _ z 2 +C 2 



£' := 



B 2 . Then for any t G R we have: 
=: a + v^T 6 = ^ |A+V ? f|+A + v^T sign^^^L^ 



and 
tg 



VW^ I _ tga + y^Tthfe _ sin(2a) + y^Tsh(25) = yzi sign^ + o^v^T 2 ) . 



1 - v^Ttg athb ch(26) + cos(2a) 
Therefore, for large \t\ we have: 



cxp 



B'(A + v^T t) - ( W 2 + P 2 ) Vx+ f Tt cotg( ^ /A+ / rTf ) + B(A+V 7\g 



cxp 



n t g (^±pn) 



(£f±il)V=I t + sign(t)(^ + 25) VA + v/^H 1+ 0(|t|- 1/2 ) 



the modulus of which is 
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e (^ + C 2 )A/2 exp 



«£±£ + 2B 2 )J lx+ ^T ltl - x + 0(1) 



Moreover 

and 

entail 



cos( v / A + V-Lt) = ch(26) cos(2a) - v 73 ! sh(26) sin(2a) 



sin( v / A + = ch(26) sin(2a) + v 73 ! sh(26) cos(2a) 



1 - cos(v / AT7 r n) - (vT+V-Tt/2) sin( v / A + v /3 T0 
= [ch(26) - cos (2a)] 2 + (a 2 + b 2 ) [ch 2 (26) - cos 2 (2a)] - 2 [ch(26) - cos(2a)] [b sh(26) + a sin(2a)] 
= [ch(26) - cos(2a) - b sh(26) - a sin(2a)] 2 + [a sh(26) - b sin(2a)] 2 

= (a 2 + 6 2 )sh 2 (26) + G(ash 2 (26)) = + O (^/\T\) 



sh 2j\\±VEph± for large \t\ . 



So far, for A < An 2 and for large |t| we have: |$«,,/3^, 2 (A + \f—\ t)\ 

= (A 2 + t 2 )3 /4 e (^ + C 2 )A/2 exp [_(«>!+£ + 2jB 2 + l )v / v^TF-A + + 0(|f|-l/2) 

= (e^V+C 2 )) |t| 3 / 2 exp [ - + 25 2 + 1 + I" 1 )) vP" 

= 0(e 2 ^ 2+ ^)exp[- (^f + 2B 2 + l/2)/p". 

Hence we can inverse the above Fourier transform for A < min{47r 2 , Ao}, and thus we 
obtain the wanted (15), which holds a posteriori for A < An 2 : 

r°° M 

e Xx qi (w,p,x,(,z)= / e-^ Itx §(\ + ^lt) — , for x > and for A < 4tt 2 . 

J-oo 27T 

Thence, for any real A < 4ir 2 , taking e < 4tt 2 — A for positive x we have 



e** qi (w,l3,x,C,z) = 0{e-") 



x 



\®(\ + £ + V-tt)\dt =0(e~ £X ), 



which entails the integrability of x (->■ e A:E gi(u>, /3, x, £, z), so that finally A > 47r 2 . o 
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